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Abstract 



Oh! 



We analyze (2+l)-dimensional Gross-Neveu model with a Thirring in- 
teraction, where a vector-vector type four-fermi interaction is on equal terms 



^ ■ with a scalar-scalar type one. The Dyson- Schwinger equation for fermion 

H ; 

self-energy function is constructed up to next-to-leading order in 1/N ex- 
pansion. We determine the critical surface which is the boundary between 
a broken phase and an unbroken one in {ac, (3c, ^c) space. It is observed 
that the critical behavior is mainly controlled by Gross-Neveu coupling 
and the region of the broken phase is separated into two parts by the line 

ac = a*(= ■^). The mass function is strongly dependent upon the flavor 
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number N for a > a*, while weakly for a < a*. For a > a*, the critical 
flavor number Nc increases as Thirring coupling /3 decreases. By driving the 
CJT effective potential, we show that the broken phase is energetically pre- 
ferred to the symmetric one. We discuss the gauge dependence of the mass 
function and the ultra-violet property of the composite operators. 
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1 Introduction 



Dynamical symmetry-breaking (DSB) plays an important role in applying 
gauge field theories to particle physics. The structure of chiral symmetry has 
been studied extensively for variety of gauge models [1 - 4] . Since the work of 
Nambu and Jona-Lasinio [1], DSB has been investigated as the mechanism 
of generating fermion masses in elementary particle physics for almost one 
generation [2-4]. Appelquist et al. [4] showed that (2-|-l)-dimensional QED 
exhibits a critical behavior as the flavor number N approachs Nc — 32 /n'^ in 
the framework of the 1 /N expansion. Such a behavior was also confirmed by 
the lattice simulation [5]. 

The search for the novel solution of Dyson- Schwinger (DS) equation 
for a fermion propagator was initiated by Johnson, Baker, and Willey [2]. 
The same equation can be derived by extremizing the effective potential 
of Cornwall, Jackiw, and Tomboulis (CJT) [6]. The CJT effective potential 
enables one to check the stabihty of non-trivial solutions of DS equation, that 
is, whether the solutions are energetically preferred to the trivial solution or 
not. 

(2-|-l)-dimensional four-fermi interaction models were shown to be renor- 
malizable in the framework of the 1 /N expansion despite its non-renormalizability 
in ordinary weak coupling expansion [7]. Such a dramatic transmutation is 
due to the fact that the composite operator acquires the large anomalous di- 
mension in strongly-correlated region. (2-|-l)-dimensional Gross-Neveu (GN) 



3 



model possesses a non-trivial ultra-violet (UV) fixed point at the leading or- 
der of 1/N expansion, which survives beyond the leading order [8]. It provides 
a ground to study DSB. (2-|-l)-dimensional Thirring model has also been 
studied in the context of the gauge structure and a dynamical mass gener- 
ation of the fermions in the 1/N expansion [9] . It is reported that Thirring 
model shows the phase transition with the strong dependence of the criti- 
cal flavor number on Thirring coupling constant, and possesses non-trivial 
UV fixed point albeit at non-perturbative order in 1/N [9]. While the gap 
equation at the leading order is local in GN model, the first non-trivial gap 
equation is non-local in Thirring model. It turns out that the gap equation is 
a relation to fine-tune the coupling as in GN model rather than the value of 
the fermion mass as in (2-|-l)-dimensional QED. The dynamically generated 
fermion mass then becomes a physical parameter which is much smaller than 
the natural cut-off of the theory. 

Thirring interaction is an another independent composite operator which 
is relevant to GN interaction in strongly coupled regime [8, 9]. Accordingly, 
it is natural for us to consider whether the above four-fermi interactions 
co-operate dynamically in the generation of fermion mass. 

In this paper, we consider (2-|-l)-dimensional Gross-Neveu model with 
a Thirring interaction, and investigate how Thirring interaction plays a part 
on the critical behavior of GN model . In Sec. 2, we introduce our model 
and derive the dressed propagators whose UV behavior improved by adding 
the gauge-fixing-like term to the Lagrangian. In Sec. 3, the DS equation is 
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constructed up to 1/N-order. We determine the critical surface which is the 
boundary between a broken phase and an unbroken one in {a, /3, N) space, 
and also present the fermion mass function in an appropriate approximation 
scheme. The critical behavior is mainly controlled by the Gross-Neveu cou- 
pling etc and the region of the broken phase is separated into two parts by the 
line etc = 0(*(— The gauge dependence of the mass function is discussed. 
In Sec. 4, we show that our nontrivial solution is energetically preferred to 
the trivial one by deriving the CJT effective potential. In Sec. 5, we discuss 
the consequences of results and the UV property of the composite operators. 



where ijj are the four-component Dirac spinors whose flavor indices are sup- 
pressed and N is the fermion flavor number. As one can see, scalar-scalar and 
vector-vector type four-fermion interactions arc parallelly introduced in the 
above equation. The 7 matrices are chosen to be anti-hermitean as follows. 



where I is 2 x 2 identity matrix , cr^'s are the Pauli matrices and j = 1, 2. 
Introducing the auxiliary scalar field a and the auxiliary vector field 



2 The Model 



Our model is given in the Euclidean version by the Lagrangian 



(1) 
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to facilitate 1/N expansion [7, 9], we can rewrite Eq. (1) as 

It is easy to see the equivalence between C and C'^'^ as follows 

£ = - In I y VaVAf, exp [- J C^H^x] | . (4) 

If we do not want to worry about what quantities should be required to be 
renormalizable, we can improve the UV behavior of all the Green's functions 
by adding to Eq. (3) the term (S^A^)^/ (2^) and postulating that the observ- 
ables are those quantities which are independent of the gauge parameter ^. 
The new theory also has a restricted gauge symmetry [9] . 

The tree propagators, after the introduction of the gauge-fixing term, 

are 

^t}{p) = h\6^,-p^p,/p^) + ^j^^{p^p,/p''), (5) 
= g\ (6) 

To find the 1 /N leading propagator, we evaluate the three kinds of one-loop 
two-point functions by the method of dimensional regularization as follows: 



= -tr J 



1 



(27r)3 '^^7 ■ k + m^" ■ {k — p) +m 

= -{Sf,u-p,Pu/p')F{p), (7) 

^/ X r d?k I 1 

= — tr / - — ; 

J [2'kY ^■k-\-m^-[k— p)+m 

= -2F(p), (8) 

d?k 1 1 



Mp) = -tr/ 



{2'kY J • k + m • {k — p) + m 
0, (9) 



where the trace is acting on a spinor space , m is the fermion mass generated 
dynamically through 1/N corrections and F(p) can be expressed without 
introducing the regulating parameter as follows [3, 4] 



F{p) = I 



m — Am^ . \ 

^+ A fZTi arcsm(^ 



(10) 



2 4vp^ V + ' 

We adapt the approximation <^ in the following analysis and in that 
case the F(p) goes to Vj?/S for the evaluation of DS gap equation in Sec. 3. 
There is no mixing between the scalar field a and the vector field A^, since 
the scalar- vector one-loop diagram vanishes (Eq. (9)). 

We can evaluate the dressed propagators by summing the chain dia- 
grams as follows: 

^'^"^^^ " 1/h^ + F{p) ^^^'^ ~ P'^P'^/P'^ + p2 +\/h'^ iPi^PJp')^ (11) 

«w = whm' '''' 

3 The Dyson-Schwinger Equation 

Wc start with the theory with a UV cut-off as was usually done in QED3+1 
[10]. By using the dressed propagators, we can construct the following DS 
gap equation [2 - 4] up to 1/N-order, 



k (27r)2 

(13) 



where 



(14) 



1-P, 



(15) 



A{p) 



T{p) 




TT- 



-1 



■27V- 
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(16) 



(17) 



E(p) is the fermion self-energy function generated dynamically and A is 
the natural cut-off of the theory. A{p) and T[p) denote the renormalization 
constants of the wave function and the vertex, respectively. They can be 
calculated perturbatively from the bare Lagrangian and their detailed 
derivations are presented in the Appendix. The second term in the right 
hand side of Eq. (13) denotes the tadpole diagram. We can not neglect the 
corrections of 1/N-order which stems from A(p) and r(p), since we want to 
consider the theory up to that order. The extra two terms denote the self- 
energy diagrams by vector field and scalar field a. The crossing diagram 
between A^,^(p) and G{p) correspond to the order of therefore, it is 

excluded from the above equations. 

Taking the trace over the gamma matrices in Eq. (13), we get the 
following DS gap equation 
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N Jo (27r3) A;2 + S2(A;) |p-/t|/8 + l//i2 

_1_ £k S(fc) 1 

~nJo (27r)3 + S2(A;) \p - fc|/4 + 



where \p — k\ = {p — ky. 

In the case of N<1 it is meaningless for us to analyze Eq. (18) since the 
1/N expansion is failed in that region. We intend to present the comparable 
feature with the previous studies [4, 5], in which it is reported that the 
fermion mass is vanished when the fermion flavor N goes to N^, = 32/71^. We 
plot the following figures in the case of N=l for illustrative purposes. 

Since the last term in Eq. (18) has the gauge parameter ^, E(p) is a 
gauge-dependent quantity. However, the fact that it is not identically zero has 
physical consequence ( i.e., chiral-symmetry breaking). Accordingly, whether 
the fermion mass is generated or not is the gauge invariant statement [9]. 

3.1 Critical surface in (ac,/3c,Nc) space 

Let us approximate T,{k) in Eq. (18) to m ft! 'E{k)\k=o- Taking the limit 
m — > 0, we can access the critical region as follows: 

{g\h^N)^iglhl,K). (19) 

In that case, the critical surface is defined by 



ac-K^ \ N^ttV iVe7r2 \ (3c J NcTT^ V a, 



c 
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+ 



^ arctan ( I , (20) 



where the dimensionless quantities are given by, 

«" = 4 ^= = 4- '''' 

Prom Eq. (20), regardless of /3c, there exists the critical flavor Nc only 
in the region of etc > a* and goes to «*(= ■^) as A^c goes to inflnity. Eq. 
(20) is plotted in Fig. 1 and Fig. 2 for two different gauge conditions, which 
are ^ = (Landau gauge) and ^ = 1 (Feynman gauge), respectively. Each 
curve in Fig. 1 and Fig. 2 corresponds to the boundary between a broken 
phase and an unbroken one. One can see that the curves of /3c = 1-0 x 10^ 
are not deformed in two figures. The other curves in Fig. 2 are shifted to 
the right for the corresponding curves in Fig. 1. 

3.2 Fermion mass function 

After taking the integrations in Eq. (18), we obtain the following equation: 

_ ^B{M) + j,C(M) + $D{M) 
l-^MM) • 

where 

M = m/A 
A{M) = 1 - Marctan(l/M) 



B{M) = 1 



7r2 



-10a -5b-— c- 26M2 + UcM^) 
3 



+ (12a - 16cM^)Marctan(l/M) 
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(25) 



C{M) = 8F{(3,M) -2F{a,M), (26) 
4 i Marctan(l/M) - aYctan{l/^) 

1 



F{a,M) = 



M2 + q;2 



— In ( ^2 J ~ arctaii(l/M) + In ( — ^ ) > 

(28) 

where we have taken the following interpolation: 

xlnx ^ ax + bx'^ + cx^, for a; G [0, 1] (29) 
with a = -2.5, 6 = 4.0, c = -1.5. 

The fermion mass functions in Landau gauge are presented for a < a* and 
q; > a* in Fig. 3 and Fig. 4, respectively. The fermion mass functions in 
Feynman gauge are presented for a < a* and a > a* in Fig. 5 and Fig. 6, 
respectively. The mass functions in Fig. 3 and Fig. 5 are strongly dependent 
on the flavor number N, while those in Fig. 4 and Fig. 6 weakly. One can 
easily confirm that the Eq. (22) is equivalent to Eq. (20) in the limit of 
0. 

4 The Vacuum Stability 

The DS equation is obtained by extremizing the CJT [5] effective potential 
with respect to fermion full propagator S{p). The CJT effective potential for 
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our model is given by 



V^'\S) = j ^tr [ln(5W -\q)S{q)) - 5^ -\q)S{q) + 1 



"^^(0) f / ^r(g)tr[%)] y 



2 ' ' V7 (27r) 

1 f d^qd^k 



2N J (27r)f 



-tr 



7^S{qh,S{q + k)A^.{k) + S{q)S{q + A;)G'(A;)J , 

(30) 



where the first term in the right hand side of Eq. (30) denotes the fermion 
one-loop contribution and the extra terms denote the two-loop contributions. 
We consider the extremal condition in the CJT effective potential as follows 



\s=s 



0. (31) 



dS{p) 

In that case, we meet the following equation: 

5(p)-^ = 5(°)-^(p) + S(p), (32) 

where, 

E{p) = VnG{0) I ^^T{kMS{k)] 

J (037.^(^)>A..(P -k)-^J ^,Sik)Gip - k). 

(33) 

One see that Eq. (32) with Eq. (33) coincides with the previous DS equation 
presented in Eq. (13). By inserting the Eq. (32) with Eq. (33) into the Eq. 
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(30), we obtain the CJT effective potential at tlie extremal propagator S{p) 
as follows 



Since the function ln(l + x) — x/{l + x) is positive for all positive x, then 
yCJT^g^ is less than or equal to zero. Hence, the energy of any nontrivial 
solutions is lower than that of the trivial (perturbative) solution E(p) = 0, 
therefore, the broken phase is always energetically preferred to the symmetric 
one. 

5 Discussion 

We studied (2+l)-dimensional Gross-Neveu model with a Thirring interac- 
tion, where a vector- vector type four-fermi interaction is on equal terms with 
a scalar-scalar type one. To solve the DS gap equation up to 1/N-order, the 
renormalization constants ( i.e., A{p), r(p)) were calculated perturbatively 
up to that order and the fermion self-energy function E(p) was approximated 
to the constant m(~ S(p = 0)). We expect that these approximations do 
not change the qualitative feature of the phase structure of our model. 

Our gap equation (Eq. (18)) contains a gauge parameter ^, and thus it 
is gauge dependent. The position A{p)'^p'^ = — E(p)^ of the pole of fermion 
propagator is a physical quantity that must be independent upon the param- 
eter. We, unfortunately, have ignored the momentum dependence of E(p) in 
our analysis, we could not present the physical mass of fermion. The numer- 




•2 + E2(p) 



(34) 



13 



ical evaluation of Eq. (18) will allow one to investigate the physical mass 
of fermion. However, as was discussed in Ref. [9], the fact that the fermion 
mass is not identically zero has physical consequence (i. e., chiral-symmetry 
breaking). Accordingly, whether the fermion mass is generated or not is a 
gauge invariant statement. The critical surfaces for ^ — 0.0 (Landau gauge) 
and ^ = 1.0 (Feynman gauge) are presented in the contour shape in Fig. 1 
and Fig.2, respectively. The curves of /3c = 1.0 x 10^ are not deformed in two 
gauge conditions, since in that case our model is equivalent to GN model in 
the bare Lagrangian level. The other curves in Fig. 2 are shifted to the right 
from the corresponding curves in Fig. 1. Such a feature is natural since it 
means that the vector bosons arc more correlated than the scalar bosons in 
the neighborhood of the phase boundary. 

The critical line etc = 0(*{— ^) in the Gross-Neveu model of the leading 
order turns into a critical surface in (0;^ /?c, -^c) space in our model (see Fig. 
1 and Fig. 2). The fermion mass function in Landau gauge are presented 
in the region of a < a* and a > a* in Fig. 3 and Fig. 4. They show that 
the fermion mass is generated dynamically in both regions. The fermion 
mass is weakly dependent on the fermion flavor number in the former, 
while strongly in the region of the latter. Such features can be reasonably 
understood. For a < a*, the critical behaviors are dominated by the tadpole 
diagram, then the theory always shows the mass generation as to all of N. 
However, for a > a*, the tadpole term does not contribute to the mass 
generation, therefore the critical behaviors are controlled by the terms of the 
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order oi 1/N and they are QED2+i-like or Thirring-like. As argued in Ref. 
[9], the effective infrared couphng weakens hke 1/N as N increase due to 
the screening effect of fermions. Therefore, the infrared couphng becomes 
so weak that the fermion condensates cannot occur in the case of > Nf.. 
All of the curves in Fig. 1 approaches to the hne a* as N increase, since 
our model in the large N limit is the GN model of the leading order (Eq. 
(18) and Eq. (20)). The fermion mass function in Landau gauge are also 
presented in the region of a < a* and o; > a* in Fig. 5 and Fig. 6. 

It is interesting to discuss the UV property [7] of composite operators ( 
i.e., {i^'^Y, and (■07^V)^) in the strongly correlated region. From the dressed 
propagators (Eq. (11) and Eq. (12)) we can see the asymptotic behavior in 
UV region (i.e., p ^ A) as follows: 

G{p) ~ -, 
P 

thus the mass dimensions of the fields in that region are [(t\uv = 1) and 
[^^ijc/y — 1- Noting the equivalence a ~ 'ipip and ~ '07iu'0 from the 
Lagrangian (Eq. (3)), we see the following facts: 

Accordingly, the renormalizability in that region is guaranteed [7]. 

Our model contains two of the most fundamental types of four-fermi 
interactions, and they co-operate dynamically in the generation of fermion 
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UV 

= 2. 



mass. In the region of a < a*, the general features of mass generation is 
not disturbed by the higher l/A'"-order corrections by the dominant tadpole 
term. In that of o; > a*, Thirring interaction plays important roles not only 
in generating the fermion mass dynamically but also in making its strong 
dependence on N. In the latter, the critical flavor number Nc is increases as 
(3 decreases. Our model, therefore, is probably a good ground for studying 
the fundamental phase transitions of the Nature. It is important to evaluate 
Eq. (18) numerically to check the stability [9] of the results. Authors are in 
process for analyzing it. 

This work is supported in part by KOSEF (Korea Science and Engi- 
neering Foundation). 
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Appendix 



We derive the renormalization constants of the wave function and the vertex 
up to 1/N-order. The former can be calculated from the self energy diagrams 
by the vector field and the scalar field, 

^' + ^= = ^ / (05>^(*)-^-S<°)(p - *) + i / (03G(*)S<°'(P - *) 
4 — 1 

(7 • p) ln(A^/p^) + finite terms. (35) 



Thus the wave function renormalization constant is given by 

'*W=l-A^'"(f)- (36) 

The latter can be calculated by the vertex diagrams by the vector field and 
the scalar field, 

1 r d^k 

V1+V2 = _y_T„A(*V„7„S(")(pi + 4)S(°'fe + *:) 

4—1 

In(AVPmax) + finite terms, (37) 



Ar3/27r2 

where p^ax largest of p{ and p| with pi and p2 the incoming and out- 

going fermion momenta, respectively. The vertex renormalization constant 
for tadpole diagram is given by 

1 3 

V{p) = ^+ lim 1n(A7pLx) 



'Kl + T^ln(AV/)). 
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Figure Captions 

Fig. 1: The projection of the critical surface on {Nc, ac) plane for various 
Pc- C is chosen to be zero (Landau gauge). 

Fig. 2: The projection of the critical surface on {Nc, etc) plane for various 
Pc- C is chosen to be one (Feynman gauge). 

Fig. 3: The fermion mass M as a function of the flavor number N for 
various a {a < a*) and P in ^ — 0. 

Fig. 4: The fermion mass M as a function of the flavor number N for 
various a {a > a*) and (3 in ^ — 0. 

Fig. 5: The fermion mass M as a function of the flavor number N for 
various a {a < a*) and P in ^ = 1. 

Fig. 6: The fermion mass M as a function of the flavor number N for 

various a {a > a*) and /3 in ^ = 1. 
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